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Abstract — Mathematics plays an important role with stereo 
data which is collected from satellite. Presently, satellite images 
can be used for various applications in civil engineering. 
Satellite stereo data plays an important role in water resources 
application by automatic generation of Digital Elevation Model 
(DEM) from aerial triangulation. 

The 3D rational function models build a correlation 
between the pixels (2D image space) and their ground locations 
(3D object space). This correlation is based on ratios of 
polynomials. 

Early in the research for the epipolar geometric 
relationship of linear pushbroom images, the exact epipolar 
curve was not easy to derive because of the complexity of the 
relationship and the difficulty in expressing it in mathematical 
form. 

Conventional aerial triangulation is reviewed. This review 
encompasses various mathematical models, self calibration 
technique, additional parameters, and the associated 
mathematical models. Mission planned satellites like IKONOS, 
QB, GeoEye and Cartosat provide the metric quality data. In 
this research work, it is proposed to use high resolution satellite 
stereo data i.e. GeoEye- 1 for creating the block setup and AT. 

Index Terms — Aerial Triangulation, Quick Bird, Digital 
Elevation Model (DEM) 

I. INTRODUCTION 

During the last two decades the space age has made a 
remarkable progress in utilizing GIS (Geographic 
Information System), Remote Sensing and 
Photogrammetry data to describe, study, monitor, and 
model the earth's surface. Improvements in sensor 
technology, especially in the spatial, spectral, radiometric 
and temporal resolution, have enabled the scientific 
community to operational the methodology. 

Early in the research for the epipolar geometric 
relationship of linear pushbroom images, the exact epipolar 
curve was not easy to derive because of the complexity of the 
relationship and the difficulty in expressing it in mathematical 
form [Otto, 1998], Subsequently, Kim [2000] derived an 
exact mathematical equation for the epipolar curve of linear 
pushbroom images using a collinearity-based sensor model.. 
However, collinearity equations can be used for linear images 
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by modifying the equations for perspective images in the 
scanning direction [Kim, 2000] . 

The end of development of remote sensing is being from 
panchromatic, multi-spectral and hyper-spectral to 
ultra-spectral with the increase in spectral resolution. On 
the other hand, present spatial resolution is impressed to 
sub-meter for example 0.5m from GEO- Eye 



> Triangulation is the process of contiguous 
densifying and extending ground control through 
computational means. 

> This operation includes establishing ground 
control points; performing interior orientation, 
Relative orientation; measuring and transferring all 
tie, check, and control points appearing on all 
photographs; and performing a least squares block 
adjustment. 

> This process ultimately provides exterior 
orientation parameters for photographs and 
three-dimensional co-ordinates for measured object 
points. 

The fundamental equation of aerial triangulation is the 
collinearity equation which states that an object point, its 
homologous image point, and the perspective centre are 
collinear 
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II. CONTROL REQUIREMENTS FOR 
PHOTOGRAMMETRIC BLOCKS 

Any block consisting of two or more overlapping 
photographs requires that it be absolutely oriented to the 
ground coordinate system. The 3D spatial similarity 
transformation with 7 parameters (3 rotations, 3 translations, 
and 1 scale) is usually employed for absolute orientation and 
requires at least 2 horizontal and 3 vertical control points. Due 
to some influences caused by transfer errors (e.g., image 
coordinate measurements of conjugate points), and 
extrapolation beyond the mapping area, the theoretical 
minimum control requirement is unrealistic. 

Controls & Block properties: Various control points 
are placed in the stereo image (GeoEye 1 ) can be seen visually 
the block properties are also with X, Y, Z. 
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Fig-4 Controls & Block properties 
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III. CONDITION AND EQUATIONS 
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Fig-2 Creation of block 



Fig-3 Block properties setup 


A. Geometric Correction Models 

In satellite imagery, a sensor model or geometric 
correction model relates object point positions (X, Y, Z) to 
their corresponding 2D image positions (x, y). Several 
sensor models can be used to correct satellite imagery: 2D 
polynomial functions, 3D polynomial functions, affine 
model, 3D rational functions, and 3D physical models 
(Tao and Hu, 2001; Fraser et al., 2002; Fraser and 
Yamakawa, 2004; Toutin, 2004a). However, for the 3D 
geometric correction of VHR satellite images, two sensor 
models are commonly used: 

(a) 3D rational functions using vendor image support data, 
and (b) 3D physical models. These two models are used in 
this work to calculate the sensor orientation for Ikonos and 
Quick Bird images respectively. Both are supported by 
OrthoEngine® from PCI Geomatica v.9.1.7 (PCI 
Geomatics, Richmond Hill, Ontario, Canada), which is the 
commercial software used here for the bundle adjustment. 

B. Collinearity condition 

The collinearity condition is illustrated in the figure 
below. The exposure station of a photograph, an object 
point and its photo image all lie along a straight line. 
Based on this condition we can develop complex 
mathematical relationships 

C. Collinearity Condition Equations 

In Fig -5 Coordinates of exposure station be XL, YL, ZL 
wrt object (ground) coordinate system XYZ 
Coordinates of object point A be XA, YA, ZA wrt ground 
coordinate system XYZ 
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Coordinates of image point a of object point A be xa, ya, 
za wrt xy photo coordinate system (of which the principal 
point o is the origin; correction compensation for it is 
applied later) Coordinates of image point a be xa’, ya’, za’ 
in a rotated image plane x’y’z’ which is parallel to the 
object coordinate system 

Transformation of (x a ’, y a ’, z a ’) to (x a , y a , z a ) is accomplished 
using rotation equations, which we derive next. 
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Fig-6 

D. Rotation Equations 
Omega rotation about x’ axis: 

New coordinates (xl, yl, zl) of a point (x’, y’, z’) after 
rotation of the original coordinate reference frame about 
the x axis by angle co are given by: 
xl — x’ 

yl = y’ cos co + z’ sin co 
zl = -y’sin co + z’ cos co 

Similarly, we obtain equations for phi rotation about y 
axis: 

x 2 = -z^in O + Xi cos O 


y2 = yi 

z 2 = Zi cos O + X! sin 0 
And equations for kappa rotation about z axis: 
x = x 2 cos k; + y 2 sin k; 
y = -x 2 sin k; + y 2 cos k; 
z = z 2 

E. Final Rotation Equations 

We substitute the equations at each stage to get the 
following: 

x = mil x’ + m 1 2 y’ + ml 3 z’ 

Wm 

i 

h 



y - m2 1 x’ + m2 2 y’ + m2 3 z’ 
z - m3 1 x’ + m32 y’ + m3 3 z’ 

Where m’s are function of rotation angles co,0 and k; 
In matrix form: X = M X’ 
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Sum of squares of the 3 direction cosines (elements of M) 
in any row or column is unity. M is orthogonal 

IV. TRIANGULATION OF GEO EYE STEREO DATA: 

In analytical photogrammetry, the collinearity condition is 
the most useful and fundamental relationship. The 
condition is expressed by the collinearity of three points: 
exposure station, object, and its photo image, in a 3D 
space. 

The image vector and ground vector are only collinear if 
one is a scalar multiple of the other [10], [11] 

Therefore the following statement can be made: 

a =KA 
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Fig 8: Geo Eye stereo data: 

In the equation, K is a scalar multiple. The image and 
ground vectors must be within the same coordinate system. 
Therefore, image vector ‘a’ is comprised of the following 
components: 

xp - xO 

a = yp- yO 
- -f - 


Where m's are the functions of the rotation angles omega, 
phi, and kappa to the object space coordinate system; and f 

is the focal length of the camera. - a are the photo 
coordinates of the image point; and xO and yO are the photo 
coordinates of the principal point that usually is known 

y y 7 

from camera calibration; J A ' A are object 

space coordinates of the object point; and 

X Y Z 

c? f - are object space coordinates of the 

exposure station. are the coefficients of the 

orthogonal transformation between the image plane 
orientation and object space orientation. The rotation 
angles omega, phi, and kappa define a sequence of three 
elementary rotations around the x, y, and z axes, 
respectively, where the xyz coordinate system is parallel to 
the XYZ object system and is constructed with its origin at 
the exposure station. 

The coplanarity condition holds that the two exposure 
stations, the object, and its image points in two image 
planes are in the epipolar plane. Because the expression of 
the coplanarity equation is difficult to calculate, it is not 
used as widely as the collinearity condition in analytical 
photogrammetry. 

V. CONCLUSION 


Where xO & y 0 represent the image coordinates of the 
principal point. Similarly the ground vector can be 
formulated as follows: 


Xp - XO 
A= Yp - YO 
Zp - zo 

Thus the following equation can be formulated as: 
a=KMA 


Theoretical and practical studies (Ackermann, 1966, 
1974 and Brown, 1979) showed that only planimetric points 
along the perimeter of the block and relatively dense chains of 
vertical points across the block are necessary to relate the 
image coordinate system to the object coordinate system. 
These measures also ensure the geometric stability of the 
block as well as control the error propagation. 

The mathematical concepts covered in this paper are: 


where, 


xp - xO 

yp ~ y° 

-f 


= KM 


' Xp - XO ■ 
Yp - YO 
. Zp - ZO . 


The above condition defines the relationship between the 
perspective center of the sensor/camera, exposure station 
and ground point. Thus the collinearity specifies that the 
exposure station, ground point and its corresponding [I] 
image point location must all lie along a straight line in a 
3-dimensional space, therefore this condition is defined by [2] 
two equations (equations 1 and 2): one equation for the x 
coordinate, and another for the y coordinate, in the image 
plane, Wolf [12]. [3] 
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> Collinearity condition equations (derivation and 

linearization) 

> Ground control for Aerial photogrammetry 

> Bundle adjustment (adjusting all photogrammetric 

measurements to ground control values in a single 
solution)- conventional and RPC based 
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